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摘 要: 本 文 研究 带 趋势 项 GARCH 误差 模型 的 参数 变 点 检验 。 在 参数 未 知 的 情况 下 ， 利 用 残 量 累积 和 
检验 以 消除 相依 序列 对 检验 结果 的 影响 :研究 表明 在 原 假 设 条 件 成 立 下 ， 残 量 累积 和 统计 量 的 极 
限 分 布 是 一 个 标准 布朗 桥 的 上 确 界 。 蒙 特 卡 罗 数 值 模拟 和 JIBM 收益 率 数据 分 析 的 结果 都 充分 说 


明了 本 文 方法 的 有 效 性 和 实用 性 。 
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1 引言 


变 点 问题 起 源 于 质量 控制 ， 从 生产 线 上 抽检 产品 以 检测 产品 质量 是 否 发 生 显著 性 波动 ， 特 
别 是 产品 质量 是 否 超 过 其 质量 检测 范围 。 当 产品 质量 超过 质量 控制 警 界线 时 ， 希 望 能 及 时 预 
警 ， 以 免 出 现 更 多 次 品 ， 这 个 质量 超过 警 界 线 的 时 刻 就 称 为 变 点 。 变 点 问题 在 工业 自动 控制 、 
金融 、 医 学 等 领域 都 有 着 实际 应 用 背景 。 对 工序 质量 的 检查 、 故 障 检测 、 分 析 汇 率 的 变化 规 
律 、 心 电 图 的 韵律 分 析 、 试 飞 数据 中 数据 跳跃 点 的 估计 等 都 可 以 归结 为 对 系统 的 变 点 研究 。 

结构 变 点 分 析 本 质 上 属于 数理 统计 和 非 线 性 时 间 序 列 分 析 范 畴 ， 它 引起 了 国内 外 众多 学 者 
的 关注 ， 并 涌现 出 大 量 关 于 结构 变 点 的 文献 ， 如 Csorgo 和 Horvath!!], Bhattcharyal?], Krish- 
naiahl3 等 。 特 别 近 十 年 来 ， 学 者 们 更 加 热衷 于 研究 各 种 包含 未 知 时 刻 结构 变 点 的 模型 。 其 
中 ， 韩 四 儿 负 ，Andrews #ll Ploberger!! 研究 了 关于 结构 变 点 的 检验 问题 ，Bail6] 则 考虑 了 关于 
结构 变 点 估计 的 一 致 性 和 收敛 速度 等 特性 。 以 上 都 是 关于 单个 变 点 的 研究 ， 而 对 于 多 个 结构 变 
点 的 检验 和 估计 的 研究 成 果 则 可 参见 文献 [7,8]。 

在 计量 经 济 学 中 ， 有 大 量 关 于 ARCH 和 GARCH 过 程 这 方面 的 文献 ， 如 : Englell, Kokos- 
zkall, Halll]. GARCH 过 程 最 主要 的 特征 是 : 序列 本 身 是 不 相关 的 ， 但 序列 每 项 平方 以 
后 的 新 序列 却 是 相依 的 。 因 为 GARCH 过程 的 条 件 方差 是 其 参数 的 函数 ， 则 可 以 通过 对 条 件 方 
差 的 检验 来 判断 其 参数 是 否 发 生 了 变化 。Andreou 和 Ghysels!!?] 提出 应 用 Inclan 和 Tiaofl3 的 
方差 平方 累积 和 统计 量 检验 GARCH(1,1) 模型 参数 的 变化 ， 而 Maekawa 等 人 [1 用 同样 的 方法 
对 ARCH 模 型 的 变 点 问题 进行 了 研究 。 尽 管 他 们 的 理论 正确 ， 但 大 量 的 数值 模拟 结果 表明 ， 
由 于 GARCH 过 程 的 条 件 方差 的 变化 会 造成 检验 的 不 稳定 ， 因 而 用 平方 累积 和 统计 量 来 检验 结 
构 变 点 会 产生 过 分 拒绝 原 假设 和 势 函 数值 过 低 的 缺陷 。 为 弥补 这 缺陷 ，Lee 等 人 15] 提出 构造 残 
量 累积 和 (RCUSUM) 统计 量 检验 GARCH(1,1) 模型 参数 的 变化 。 这 是 因为 基于 残 量 的 检验 能 
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够 消除 序列 彼此 相依 的 影响 ， 从 而 改进 检验 的 效果 。 但 是 由 于 一 般 的 GARCH(p,9) 模型 具有 
更 复杂 的 概率 结构 ， 所 以 把 文 [15] 的 结果 直接 推广 到 GARCH(p, q) 模型 场合 存在 一 定 的 难度 。 
本 文 应 用 Berkes, Horvath 和 Kokoszkal16| 对 GARCH(p, 9) 模型 研究 的 一 些 结果 ， 对 带 趋势 项 
的 GARCH 误差 模型 的 参数 变化 进行 残 量 累积 和 检验 。 
2 ”模型 和 假设 

本 节 沿 用 Berkes 等 人 86] 的 定义 和 假设 ， 考 虑 模型 


Ye 三 4 十 十 Ti Me = heét, 


Pp q 
he = wt > amit Y bhy, (1) 


i=1 j=1 
其 中 mh 是 一 个 GARCH(p, gq) 过 程 。 定 义 参数 向 量 
T = (u, p,0) = (4,8,w,0,B1), 1" = (u*, 8", 6") = (u*, B*,w*, a, Bis), 
Hrn JRA ER A ER St HE SCF 
Ho: BR nr = (jp,B,w, Qa1,B1,… ) 保持 不 变 ， 
Ay: 参数 在 时 刻 k* Mx = (u,b, w, a1, Bi,---) 变化 到 zr* = (u*, B*,w", at, Bi). 
ENRI c; (—00 < t < 00) 是 独立 同 分 布 序列 ， 令 


È = 1 T 1 loz ĝ ne 
T(u) = Troal og (u) + ee 


其 中 : 
be(u) = co(u) + 2 ci(u)772 i, 


i=l 
图 数 c;(u) (0 < i < co) 在 文献 [16] 中 定义 如 下 。 对 于 w= (z,51, ,spyt1,… te) Æq > p， 
则 
co(u) = 7/(1— (ti +--+ +tg)), 
ci(u) = s1, 


c2(u) = s2 + tıcı (u), 


Cp(U) = Sp + t1Cp-1 (u) + tocp_a(u) +--+ + tp-161 (u), 


Cp+1(u) = ticp(u) + tocp_i(u) +++ tpc1 (u), 


calu) = tyeg—1(u) + tocg—2(u) +--+ +tg—101(u). 
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若 g <p， 则 ci(w) (0 <i< co) 定义 为 
co(u) = 2/(1— (tı ++ to)), 
ci(u) = $1, 


co(u) = s2 + tici (u), 


Calu) = Sq + t1Cg_1(u) + t2cq-2(U) +--+ + ta-1061 (0), 


co+Hli(u) = sg+1 + tica (u) + tocg—1(u) +++ + taci (u), 


Cp(u) = Sp + t1Cp—1(u) + tacp_a(u) + +++ + tgCp-q (u). 
# i> R=max(p,q), H 
ci(u) = tici_1(u) 十 toci_2(u) +--+ tgCi—q(u). 


ARIE GARCH(p, q) 过 程 具 有 唯一 平稳 解 ， 引 入 以 下 符号 。 令 


TT = (Ai t ayer, Ba, aaa: | Bq-1) € Rv}, 
sr = (Ef, ,0) € RT, 
a = (az, ,ap 1) € RP. 


不 失 一 般 性 ， 假 设 min(p,g) > 2. LM (p+q—1) x (p+ q—1) HARE Ar: 
Tr ba Q ap 
下 
OO. “Bek 0 


FOP -1 M Ip- aN g 一 1 维 与 9 一 2 维 的 单位 矩阵 。d x d 维 矩阵 M 的 范 数 定义 为 


|M|| = sup{||@X'||a/l|Xla : X € RY, X #0}, 
其 中 :是 Re 空间 的 Euclid 范 数 。 

序列 {AT, -oo < T < co} 的 最 大 Lyapunov 指数 定义 为 
YL = inf 元 Blog | Ao-A1 sees Ar|l. 


1<T<oo 
设 
E(log ||4ol) < 00, (2) 
Bougerol 与 Picardl7,19 证 明 若 上 式 成 立 ， 则 GARCH(p,9) 过程 有 唯一 平稳 解 的 充分 必要 条 件 
是 
YL < 0, (3) 
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注意 到 这 个 条 件 隐 含 着 D +--+ + Bq <1。 
现在 介绍 参数 空间 的 一 些 条 件 。 令 0 < u< Us 0 < Po K 1, qu < Po» 定义 


U = {u: ti +ti +- +t < po, U< min(2, 81, 82,-++ ,Sp,t1,.** , tq) 


< Imax(Z; 851, 82,* Sp b1, sta) <u}. 


假设 9 e U， 此 条 件 意味 着 U 中 的 每 一 个 元 素 可 以 作为 GARCH 过 程 的 参数 。 
为 使 参数 9 是 唯一 确定 的 ， 需 以 下 两 个 条 件 : 


(i)” 实 系数 多 项 式 aiz +a? +--+ +apz? 与 1 一 Diz 一 Box? 一 .… 一 Box? HR, 
(ii) sl 是 非 退 化 的 随机 变量 。 
为 估计 wo(w)/wo(90) KIHE, Berkes AS) BEEE k > 0 使 得 


lim, t-“P(e? < t) =0. (4) 
对 sl 还 需 附 加 以 下 条 件 
Ee? = 1， (5) 
存在 5 > 0， 使 得 
Ele,|*** < 00. (6) 
Berkes 等 人 06] 证 明 若 不 存在 变 点 条 件 ， 且 条 件 (2)-(6) 成 立 ， 则 TW2(6 — 0) HSCF IE 
态 分 布 ， 并 且 有 以 下 两 个 引 理 成 立 。 
5 [38 2.1106) MER u = (T, $1,, 82, Ba 1 Sp, b1; nae ita), u* = (x*, s1,, 83, een , Sp tt, os sta) € 
U, & 


* 


ija Glu) o ; 
Ciut < ci(u), cilu) < Copy", alu) = Co( max 3 V 1) ， 0<i<o. 


其 中 01, Co, C3 为 正常 数 。 
318 2.209 设 0e U， 则 有 


Jo e ðco(u) , 
; = <j< 
| Ox Evers 0s; 0, 1l<j<p; 
colu) u ‘ Ic, (u)| | | 
<=. l<j<q: A fh . 
at; Beers SISA aG) <C4i, 1<i<oo 
其 中 C4 是 常数 ， 
,~ /Oci(v) Oci(u) Oci(u) Oci(u) Beci(u) 
A Bae a a 
3 ”主要 结果 


本 节 给 出 GARCH(p, q) 模型 参数 变化 的 残 量 累积 和 (RCUSUM) fH. Ky, t= 1,2,.… ,7 了 T， 
是 样本 数据 ， 参 数 9 的 估计 定义 为 


6= arg max{Ly(u) :weU}. 
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条 件 方差 
= co(0) + X c (0)h?_, 
j=l 
的 估计 为 
= co(ĝ) + DEOL ; 
j=l 
其 中 和 是正 整数 序列 ， 满 足 条 件 : 当 了 一 co 时 ， 和 一 0 HA?/T — 0. 
则 残 量 e 的 估计 为 


入 
= (6) + > (6) A2_5, th = ye — Ê — Bt, és = the/he, 
j=l 
其 中 有 一 = O,(T-"/2) F Â — p = O(T-3/2)。 基 于 et 的 累积 和 统计 量 定义 为 
ne. SE ee 
mele pL 
Hc, 是 有 界 平稳 的 随机 变量 ， 且 7? = Var(c?), FREE 的 累积 和 统计 量 定 义 为 
1 i ire 
fee ames, | oa - zal 
这 里 
a4 2 a2 2 
= as (= 2 ae 
定理 3.1 若 引 理 2.1 和 引 理 2.2 成 立 ， 对 任意 的 0 < po < 1， 存 在 一 个 正 整 序列 和 = Ap 


满足 当 T 一 oo 时 ， 有 和 一 oo0， 和 /VT 一 0 条》_ pi? < oo， 且 序列 {(m,hi)} 是 强 混合 序 
列 ， 则 当 原 假设 Ho 条 件 成 立时 ， 有 


Ê, & sup |w? (v)], 
0<v<1 


HP v= k/T, W(v) X [0,1] 上 标准 的 布朗 桥 。 


证 明 注意 到 
a2 a Ay\2 

2 ht 

awa at t+2 

t h2 he h? + gz Pui - b) m(Ê — n) 

+AA- n)(B— B)e+ (i-u? + PB —p)] = 4 EY 

=3 

2 2 _ h2)e2 2 _ 7,2)2,2 2 

m _ Fr a -hijet , (he —he) et _ 2 

he ae RT ARP =e + D Vat- 


要 证 明定 理 3.1， 实 际 上 等 价 于 证 明 下 面 结论 的 成 立 。 


Pa = Fp | So V- FD Vel = ol, 1 一 1 7. (7) 
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第 26 郑 
我 们 首先 处 理 Vi,:。 注 意 到 
oo 入 
he — he = (co(8) — co(6)) + > cOn +Y (cr(b) cÂ; 
J=A+1 j=l 
入 4 
十 > Ân; 一 名 让 = YO Wis 
j=1 i=1 
根据 强 混合 序列 的 不 变性 原理 ， 有 
1 k ye? E2 k we E? 
i bse Do (a Er) ~ 7D (Ga Erg) = O0. (8) 
由 引 理 2.1， 引 理 2.2 及 中 值 定理 ， 有 
1 A ke2 ke? 
ae 
= |(6 - 6) TENDE Silá Co)| = op(1) : Op(1) = op(1). (9) 
同样 ， 根 据 强 混合 序列 的 不 变性 原理 ， 有 
1 Erm; k Es 
Tp Sth — hy St) — og 
所 以 
1 Pet ke em 
CC ao T h? *)| 
sol SD 0)) = Op (of *P/4) = op(1). (10) 
了 一 入 十 1 
同 理 可 证 ， Wt 满足 下 式 
1 e?ne k wn co 
= 046) (> 2-2 j 
Tm (6) sD h? T h? )| 
入 ` A 入 
= (X — DIGI) = Op( 9-10 — Alile) 
j=l j=l 
‘ j/ AS 
q\ _ q\ _ 
SC Op ) = Oo dm )= an 
最 后 ， 只 需 验证 Wat 满足 
1 We? k Wate = 
i Bs | 应 T2 ig |= ar) 
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注意 到 


c;(ô) [m2 ; 一 Ñe] = c;(6)[2 (A — uym- + 2(8 — B)tne-; 


—2(fi — „(Ê — 8)t — (à - p}? — (Ê — 8)? ] -Pa 


为 证 (12)， 需 验证 


1 amen QitE k È QitE 
= itt a,tE¢ 
到 = 1 =] 
Biel Oa PER ae [eG isbot 
因为 
1 ?1 一 7E ke MjE 
max t-jSt t—jEt 
| 2 SSS 2 = O,(1), 
VT 1<k<T tel h T ram hi 


FEA- u= O0(T"?), Ĝĝ- p= OT), 所 以 


HE Lowa (DT L Ey Ho a) 


= o (1Y) = M + op(1) = 0p(1). 


j=1 


又 因为 


k T 
1 tm—je? k tme—j€? 
da: at = 0,(7). 
T (> h2 T h? ) p 


t=1 t t=1 


同 理 可 得 


1 mx > WB- A ky sms)| 
VT 1SkST 和 有 T R 


入 入 
= (TE ô) = (T25 cs(0)) + op(1) = 0p0). 
j=l 


j=l 
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(13) 


对 于 Qat, Qat, Q5t KE WARE Qe Qoe 的 证 明 是 一 致 的 ， 故 略 。 这 样 我 们 证 明 


下 Rit = = op(1)。 


下 面 我 们 来 证 明 Roe = op(1)。 注 意 到 处 > 0， 那 么 我 们 总 可 以 找到 一 个 正常 数 M > 0， 使 


得 如 > M。 因 为 


4 
Ta 》 Wit, 
i=1 
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那么 先 证 明 W2, 满足 下 起 


1 g 27 mje k T Me-—jEt 
ea )) ep hehe -TÈ heh? )| 
1 2 me?2 k Ane 
= VT 1cker Laz amy, ( 2 a0) (> T - TŽ EE e) 


= 0,(( D o0) ) =0,(( 3 À) = op(1). (14) 
j= 和 +1 j= 和 +1 


同 理 可 证 W, 满足 下 式 


入 k 4 22 T V4 22 
万 [60-60) Le TA eee | 
入 k 
< ae ee Fe max 1/7 (> (ci(9) 加 cj(O) (过 we ; > te) 
= 0,((\0- 45666)’ = 0,( (6 - AY seul J 
j=1 j=1 
0,( (Ie - Dy ÀN) =0 o (ŽEN) = 0p(1). (15) 
j=1 
由 于 
=a XY whe? SA < Ww rE? 
t=1 t=1 
要 验证 Ro = 0,(1)， 还 需 验证 
FE W2? = op(1)，  i=1,4. (16) 


由 引 理 2.1， 所 以 co(9) — co(ô) = (8 — 8) - colo) = Op(T-1/2)， 并 有 


T 入 
oS (Sale )(ne ~ 9). Ef = 0p(1). (17) 


1 j=l 


Sl- 
IMs 
N 
V 
i 
引 - 

Il 


前 面 已 证 
5 
ci (Â) (n; = fez) = > Qie 
i=1 
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则 只 需 证 明 下 式 成 立 
1 工 24 
VE x (X Qie) Et = 0,(1), i=1,... ,5. (18) 


t=1 j=l 


注意 到 


元 (Es Oym-s(—n)) e? 


j= 
=0 (LY ro(So@)’) =0 exes (0) ) + 1) = 0,(1 19 
= OTE a =r JT atid Op(1) = op(1). (19) 
E 容易 证 明 


F Ds (X ôm- - yt)” e? 


t=1 j=l 


— 


T 


5 0,( 元 Ser 2«0)) = (zl E0’) + op(1) = op(1). (20) 


同 理 ， 对 于 Ose, Qai @s 的 证 明 本 质 上 与 Qi Qoi 的 证 明 也 是 完全 一 致 的 ， 我 们 再 次 省 
略 其 证 明 步 又。 这 样 我 们 完成 了 Roy = op(1) 的 证 明 。 
类 似 于 (16) 的 证 明 ， BERRIES ia 有 


T Lyw Her =op(1),  i=2,3. (21) 
联 立 (16) #1 (21), MI 
T 
= DOA — be)2e? = of (1). (22) 
t=1 


我 们 还 要 证 明 当 i = 3,… ,7 时 ， 有 Ri = op(1)。 很 容易 证 明 Ret = op(1)。 此 外 ， 我 们 完 
全 可 以 用 证 明 Rit = op(1 ) 和 Ras = = S ATER EE a = op(1), i = 3,4,5,7. ERA 
就 证 明了 (7) 式 ， 也 意味 着 


k T k k T 
E apd a= p my De- ,| + op(0). (23) 
由 泛 函 中 心 极限 定理 
k 1 Ño 2) d 1 2 2, d 
万 Vir Ye — Eez) — vW (1), JT- DG 一 五 cf) — W (v), 
t=1 t=1 
成 立 ， 再 由 连续 映射 定理 ， 则 


1 k < 
a 2 2 
| | 


k T 
__1 2 n 天 2 z| < 0 
= Tay MAY, | Lie as eae 2e -~ Ee?)| -+ pee, IW |. 24) 
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要 证 明定 理 ， 最 后 我 们 还 要 证 明 方差 估计 的 一 致 性 : 72 5 7? = Var(e?) = Eet 一 
注意 到 、 
oI -n 
根据 前 面 证 明 ， 我 们 知道 ES 满足 


因此 ， 由 (22)，(25) 和 (26) 式 ， 有 


t=1 t=1 
这 样 就 证 明了 

1 es P 

T t=1 z Hel 


同 理 ， 根 据 (22), (25) 和 (26) R, A 


T 
1 人 
$3 ê = e?) 2g TM > Ate = 82)2c4 + op(1) 


1 
S (Bs ed): (7 D0 — he) e?) + op( ) = 0p(1). 
从 而 有 
Ipea 二 ise gis a 
TE t t T E t t 了 全 
2 £ SF 
oe dé =e) Fa Le = O,(1) 
因此 
工 T 
i 
Dreo < (2 Da a-a "(Lye sey)” 
ae t=1 t=1 
= op(1) + Op(1) = op(1). 
这 样 就 证 明了 


DE -Zs Bet. 


联 立 (30) #1 (31) xt, 72 5 72。 再 鉴于 (23) 式 和 (24) 式 ， 则 定理 得 证 。 
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(25) 


(26) 


(27) 


(28) 


(29) 


(30) 


(31) 


第 6 期 金浩 等 : 带 趋势 项 GARCH 误差 模型 参数 变化 的 残 量 检验 1079 


4 ”数值 模拟 和 实例 分 析 


本 节 首 先 用 数值 模拟 说 明 用 残 量 累积 和 方法 检验 带 趋势 项 GARCH 误差 模型 参数 变化 的 可 
行 性 。 为 了 更 好 的 分 析 检验 统计 量 Ry 的 性 质 ， 我 们 考虑 GARCH(1,2) 模型 


Yọ = pt Pt+m, 
m = hier, 
h? = w+ ana + Bihi + Boh?_». 


我 们 用 蒙特 卡 罗 方 法 研究 在 大 样本 容量 条 件 下 的 残 量 检验 统计 性 质 。 令 st 是 独立 同 分 布 
的 正 态 随 机 变量 序列 。 设 置 变 点 时 刻 为 kr* = T/2， 即 参数 从 r = (yu, B,w,a1, Ai, b2) Bb 
到 天 = (pj*,B*,w*,at, Br, 访 )。 样 本 容量 分 别 依次 取 为 了 = 500, 1000, 2000， 且 检验 的 显著 性 
水 平 为 0.05。 此 外 取 入 = [logTJ?， 并 用 1000 次 试验 中 拒绝 原 假设 的 百分数 作为 检验 水 平 值 或 
经 验 势 函数 值 ， 模 拟 结果 见 表 1- 表 3。 


Al: ”经验 水 平 值 和 经 验 势 函数 值 ， 其 中 7 = (0.02, 0.28, 0.2, 0.2, 0.3, 0.3) 


x n = 500 n = 1000 n = 2000 
(0.02, 0.28, 0.2, 0.2, 0.3, 0.3) 0.033 0.036 0.043 
(0.02, 0.28, 0.3, 0.1, 0.1, 0.4) 0.481 0.742 0.907 
(0.02, 0.28, 0.1, 0.3, 0.5, 0.2) 0.557 0.843 0.942 


表 2: ”经 验 水 平 值 和 经 验 势 函数 值 ， 其 中 7 = (0.02, 0.28, 0.2, 0.2, 0.3, 0.3) 


T* n = 500 n = 1000 n = 2000 
(0.02, 0.28, 0.2, 0.2, 0.3, 0.3) 0.036 0.038 0.049 
(0.03, 0.38, 0.2, 0.1, 0.1, 0.5) 0.494 0.754 0.917 
(0.04, 0.48, 0.2, 0.3, 0.5, 0.1) 0.596 0.859 0.958 


表 3: 经验 水 平 值 和 经 验 势 函数 值 ， 其 中 六 = (0.02, 0.28, 0.2, 0.2, 0.3, 0.3) 


a n = 500 n = 1000 n = 2000 
(0.02, 0.28, 0.2, 0.2, 0.3, 0.3) 0.041 0.048 0.060 
(0.03, 0.38, 0.3, 0.1, 0.2, 0.5) 0.681 0.892 0.957 
(0.04, 0.48, 0.1, 0.3, 0.4, 0.1) 0.657 0.903 0.982 


从 3 个 表 中 可 以 看 出 : 在 原 假 设 成 立 的 情况 下 ， 残 量 累积 和 (RCUSUMD) 检验 统计 量 没有 出 
现 过 分 拒绝 原 假 设 的 现象 ， 同 时 随 着 样本 容量 的 增加 ， 经 验 水 平 值 和 经 验 势 函数 值 都 有 增加 。 
在 不 同 的 样本 容量 下 ， 经 验 水 平 值 也 都 较 接 近 显 著 水 平 0.05， 即 检验 方法 在 原 假设 成 立时 犯 第 
一 类 错误 的 概率 都 很 小 。 反 之 ， 如 果 当 备 择 假设 成 立时 ， 出 现 了 一 些 有 趣 的 现象 。 我 们 发 现 
如 果 参 数 0 = {w, a1, Bi1, 62} 变 化 的 幅度 越 大 ， 则 经 验 势 函数 值 越 大 ， 说 明 样 本 中 出 现 变 点 的 
概率 也 越 大 。 虽 然 其 他 参数 jy 和 6 对 经 验 势 函数 值 的 影响 相对 于 参数 0 = {w, a, b, Go} 594, 
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但 这 并 不 能 否定 这 些 参数 对 变 点 检验 的 影响 。 究 其 原因 ， 显 然 这 是 因为 样本 条 件 方差 主要 是 
由 GARCH 的 参数 所 决定 。 总 之 ， 通 过 数值 模拟 的 分 析 充 分 说 明了 用 残 量 累积 和 统计 量 检验 带 
趋势 项 GARCH 误差 模型 的 参数 变化 的 合理 性 。 

为 了 进一步 说 明 残 量 累积 和 统计 量 的 特性 ， 我 们 用 IBM 从 1961 年 5 月 17 号 到 1962 年 7 月 2 
号 的 收益 率 数据 做 实例 模拟 。 数 据 的 原始 序列 见 图 1， 显 然 % 是 非 平 稳 的 。 而 对 数 差 分 的 序 
Fn ye 一 Iny_1 见 图 2， 并 对 此 序列 应 用 平稳 性 游程 检验 ， 得 到 游程 数 + = 239， 其 均值 Er = 
182.42， 方 差 D(r) = 7960.86， 在 检验 水 平 0.05 下 ,得 2 = (r—Er)/yD(r) = 0.6365 < 1.96, 
故 认为 对 数 差分 序列 是 平稳 的 。 尽 管 如 此 ， 我 们 从 图 2 中 发 现 对 数 差分 序列 的 波动 是 不 均匀 
的 ， 我 们 猜测 这 可 能 是 原 序列 的 分 段 平稳 造成 的 。 在 检验 水 平 0.01 下 ， 用 残 量 累积 和 统计 量 对 
原始 序列 检验 ， 接 受 变 点 存在 的 假设 并 得 到 变 点 估计 天 = 234， 其 中 变 点 估计 的 统计 量 为 


k= min {k: [Rel = max |R}, 


并 用 BIC 准则 确定 GARCH 过 程 的 阶 数 ， 结 果 见 表 4。 


#4: 用 GARCH(1,2) 模型 拟 合 IBM 股 票 价格 原始 数据 的 结果 


H B w a Bı Ba 
[1,234]  0.07(0.356) 0.23(0.671) 17.89(2.81) 0.36(0.19) 0.38(0.06) 0.04(0.08) 
[235,350] 0.46(0.443) -0.66(0.713) 49.22(7.22) 0.44(0.21) 0.45(0.09) 0.07(0.13) 


1 ro 1 1 
o sO 100 150 200 250 300 


图 1: 数据 原始 序列 图 


本 a n à SEE 1. 
o 50 100 150 200 250 300 


图 2: 数据 对 数 差分 图 
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在 变 点 时 刻 之 前 ， 即 从 时 刻 上 = 1 Blt = 234 的 模型 为 
Y; = 0.07 + 0.23t + m, 
nt = heét, 


h? = 17.89 + 0.3672_, + 0.38h2_, + 0.04h2_,. 


而 在 变 点 时 刻 之 后 ， 即 从 时 刻 上 = 235 到 t = 280 的 模型 为 


5 


Y, = 0.46 — 0.66t + t, 

M = hires, 

h? = 49.22 + 0.44n2_, + 0.45h?_, 十 0.07h?_». 
结论 


本 文 沿 用 由 Lee 等 人 (1 引 提 出 的 残 量 累积 和 统计 量 检 验 带 趋势 项 的 GARCH 误差 模型 的 参数 


变化 。 研 究 表明 如 果 原 假设 条 件 成 立 ， 残 量 累积 和 统计 量 的 渐 近 极限 分 布 是 一 个 标准 布朗 桥 的 
上 确 界 。 证 明 过 程 中 ， 我 们 用 到 了 由 Pligradl19 与 Carrasco 和 Chenl?0 证 明 的 关于 强 混合 序列 
的 不 变性 原理 。 数 值 模拟 结果 表明 ， 残 量 累积 和 检验 的 经 验 水 平 值 和 经 验 势 函数 值 都 较 接 近 预 
先 设 定 的 显著 水 平 。 总 之 ， 对 于 GARCH 模型 参数 变化 的 检验 问题 ， 残 量 累积 和 统计 量 不 失 为 
一 个 有 效 的 方法 。 
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Residual CUSUM Tests for Parameters’ Change in GARCH Error 
Models with Deterministic Trend 
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Abstract: This paper analyzes the test problem for parameters’ change in GARCH error models 
with deterministic trend. The residual cusum test is proposed when the parameters are unknown, 
which can weaken the influence of dependent observations. It is shown that the asymptotically limiting 
distribution of the residual cusum test statistic is the supermum of a standard Brownian bridge under 
the null hypothesis. In order to verify our conclusion, we carry out a Monte Carlo simulation and 
examine the return of IBM data. The results from both simulation and real data analysis support our 
claim. 


Keywords: residual cusum test; invariance principle; Brownian bridge; change point 


